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PREFACE 


This atlas presents atmospheric general ;irculation statistics derived from twelve years 
(1979-1990) of daily National Meteorological Center (NMC) operational geopotential 
height analyses; it is an update of a prior atlas using data over 1979-1986. These global 
analyses are available on pressure levels covering 1000-1 mb (approximately 0-50 km). 
The geopotential grids are a combined product of the Climate Analysis Center (which 
produces analyses over 70-1 mb) and operational NMC analyses (over 1000-100 mb). 
Balanced horizontal winds and hydrostatic temperatures are derived from the geopotential 
fields. This atlas presents monthly means, along with daily and inierannual variabilities, 
of various circulation statistics derived from these data: zonal mean temperatures ar 4 
winds, poleward heat and momentum fluxes, Eliassen-Palm flux divergences, and 
stationary and transient planetary wave statistics. Meridional cross sections are presented 
for each month, along with Northern Hemisphere and Southern Hemisphere polar 
projections at various pressure levels for January, April, July, and October. 
Climatological latitude-time and height-time cross sections of zonal mean winds, 
temperatures and wave fluxes are presented to define seasonal variability, and actual time 
series for each year are included to show true daily variability. Finally, tables of monthly 
average zonal mean winds and temperatures are given for reference. This atlas is 
presented for a concise climatology of dynamical quantities throughout the troposphere 
and stratosphere, as a climatology of planetary wave statistics, and to alert potential users 
to the availability of these data. 


Acknowledgments: This updated version of this atlas was completed in the Atmospheric 
Chemistry Division at NCAR, with support from NASA Grant W-16215. The majority 
of the work was actually completed in connection with the prior version, under support 
from the Advanced Study Program at NCAR. Constructive reviews were supplied by 
Rolando Garcia and Kevin Trenberth. Marilena Stone and Urusla Rosner helped with 
preparation of the manuscript. NCAR is sponsored by the National Science Foundation. 
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1. DATA DESCRIPTION 


ORIGIN 

The data used here are horizontal wind and temperature fiet '? derived from daily 
geopotential height grids covering 1000*1 mb for the years 1979-1990. Grids at pressure 
levels 1000, 8S0, 700, 500, 400, 300, 250. 200, 150, and 100 mb are 1200 GMT output of 
the Global Data Assimilation System (GDAS) used at the National Meteorological Center 
(NMC). This system is described by McPherson et al. (1979). These analyses result 
from optimal interpohti jn of all available ob.>ervations with first-guess fields provided 
by the NMC prediction model. 

Stratospheric grids at 70, 50, 30, 10, 5, 2, 1 and 0.4 mb are 1200 GMT daily 
operational analyses produced by the Climate Analysis Center (CAC) section of NMC. 
These data are produced using available radiosonde observations (RAOBS) in a modified 
successive correction (Cressman) analysis, with satellite-derived tempera. lire retrievals 
used as first-guess fields (see Gelman et al., 1986, and references therein). The satellite- 
derived temperature retrievals are obtained principally from the TIROS Operational 
Vertical Sounder (TOVS) series, which have been operational since late 1978. TOVS is 
comprised of three radiometers which observe the stratosphere and troposphere (Smith et 
al., 1986): the High Resolution Infrared Radiation Sounder (HIRS-2), the Microwave 
Sounding Unit (MSU) and the Stratospheric Sounding Unit (SSU). Because the 
uppermost analysis level at 0.4 mb lies substantially above the peak of the uppermost 
satellite radiometer weighting function (SSU channel 3, which peaks near 1.5 mb), data at 
this level are of suspect quality and not used in this atlas. As a note, CAC produces 
separate analyses of temperature and geopotential height at stratospheric levels. The 
temperatures used in this atlas are not based on their direct temperature analyses, but 
calculated hydrostatically from the height fields. 
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Other recent compilations of stratospheric climatological data include Hamilton 
(1982), Geller et al. (1983) and Wu et al. (1987), v ho use subsets of the NMC/CAC data 
used here, and geostrophically-derived winds. Labitzke et al. (1985) provide a 
climatology of monthly mean zonal mean temperatures, zonal winds and stationary 
planetary wave amplitudes and phases up to 80 km (0.01 mb). Their analyses are based 
on two years of Nimbus 5 Selective Chopper Radiometer (SCR) combined with three 
years of Nimbus S Pressure Modulated Radiometer (PMR) inverted radiances. Marks 
(1989) presents balanced winds and stationary wave fluxes derived from these data. 
Fleming et al. (1988) extended this climatology up to 120 km (0.0001 mb), providing 
useful cross sections and other derived quantities. Comparisons between these data and 
NMC/CAC derived statistics show reasonable agreement; some differences are found 
between calculations based on geostrophic versus balanced winds. This atlas furthermore 
includes transient wave: statistics, and documents the daily and interannual variabilities of 
dynamical quantities over 1000-1 mb. Apait from the changes in satellite retrievals and 
analysis schemes described in the next section, it is clear that twelve years of data is 
insufficient to quantify the interannual variability of the stratosphere (Labitzke and 
Naujokat, 1983, for example, examine data from 1955-1982); the objective here is to give 
qualitative and structural information. 

CHANGES IN 0£ERATIQ N ALj\NAL Y S1S SYSTEMS 

Because statistics of interannual variability are presented here, it is important to 
be aware of changes in the analysis schemes. Both the NMC 1000-100 mb analyses and 
the CAC 70-1 mb data are operational products, and as such are continually undergoing 
changes in both processing and sources of data in order to produce the best product. 
Changes in the NMC GDAS assimilation system and resulting discontinuities in products 
over 1979-1986 have been summarized in Vrenberth and Olson (1987). Important 
modifications in the NMC forecast model used to generate the first guess fields include a 
change from gridpoint to spectral model and the use of nonlinear normal mode 
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initialization (May 1980), various changes in horizontal and vertical resolution, and a 
host of changes in data use criteria. Trenberth and Olson's analyses suggest lhat these 
changes have the largest impact in the tropics and for fields such as relative humidity; the 
impact on extratropical geopotential height analyses are probably smaller. One region 
where there were obvious problems in the NMC analyses was in the troposphere over 
Antarctica prior to May 1986 (Bonner et al. , 1986); correction of this problem leads to 
large interannual temperature variability that is purely an artifact of the change in 
analyses. 

Changes in the stratospheric 70-1 mb analyses include modifications of the 
satellite-derived first guess fields, and differing criteria for the inclusion of RAOB data 
over 70-10 mb; these are summarized in Gelman et al. (1986). Direct inclusion of RAOB 
daia in the SH analyses was discontinued after 17 October 1980. Satellite temperature 
retrievals prior to 17 October 1980 were produced via linear regression of Vertical 
Temperature Profile Radiometer (VTPR) and SSU channels 1 and 2 measurements. 
From 17 October 1980 to the present, temperature retrievals from TOVS have been used. 
A series of TOVS instruments have been put into orbit aboard a succession of operational 
satellites; these instruments do not yield identical radiance measurements for a variety of 
reasons, and derived temperatures may change substantially when a new instrument is 
introduced (Johnson et al ., 1986; Nash and Forrester, 1986; see also Chapter 6 of WMO, 
1988). Gelman et al. (1986) have compared the operationally derived temperatures with 
co-located radiosonde and rocketsonde temperatures, and find systematic biases for each 
different satellite system of order +/- 3-6 K in the 3-0.4 mb analyses. They suggest that 
corrections be applied to the satellite-derived analyses at and above 5 mb to remove these 
biases. Such corrections have not been applied to the data used here, and a component 
of the calculated interannual variability in temperatures at these levels is due tn changes 
in data r equisition. This can be clearly seen, for example, as localized maxima near S mb 
in zonal mean temperature interannual variability maps. Note, however, that such biases 
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are not strong functions of horizontal position, and derived horizontal winds are probably 
not as strongly biased by the introduction of different satellite instruments. Therefore, 
interannual variability of the upper stratosphere in these data is probably best viewed in 
terms of zonal wind variability. 

The stratospheric CAC analyses are produced separately for each hemisphere on a 
65 x 65 polar stereographic grid. Although all the original CAC 65 x 65 stratospheric 
grids contain full fields, data equatorward of 20* in the Northern Hemisphere was purely 
extrapolated from higher latitudes until 29 February 1984. Following this time, the first 
guess fields have been archived over 0-20* N. To provide complete analyses here, the 
geopotential fields prior to 1 March 1984 have been linearly interpolated over 0-20* at 
and above 70 mb. 

GEOPOTENTIAL GRID PROCESSING 

The geopotential height data over 1000-1 mb have been archived via zonal 
Fouriei components, up to zonal wave number 12, on a 40-point Gaussian grid in latitude 
(with approximate spacing of 4.5 degrees). This latitude resolution is that associated v.,th 
rhomboidal spherical harmonic truncation at wave number 15 (R15). The stratospheric 
data are obtained from NMC on 65 x 65-point polar stereographic grids; these were 
interpolated via a Bessel 16-point two dimensional scheme to a latitude-longitude grid, 
followed by zonal Fourier analysis. Some of die 1000-100 mb NMC data were obtained 
from the 2.5° x 2.5° latitude-longitude NMC grids archived at NCAR. 

Missing data (less that 5% of the grids) have been linearly interpolated in time, 
and some obviously bad data have been replaced by interpolated values. The global data 
are stoned by month %itd year on the Mass Storage System at NCAR. Prospective users of 
the data are welcome, and may simply contact the author. 
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TEMPERATURE. WIND AND WAVE FLUX DERIVATIONS 


Temperatures are derived grom the geopotential grids via the hydrostatic relation: 


T = 



( 1 ) 


Here d> is the geopotential height, In p is the logarithm of pressure, and R the gas constant 
(287 $-2 K' 1 ) Although simply called temperature here, this expression actually 

results in the virtual temperature (the temperature of dry air having the same density and 
pressure as the true, moist air). 


Zonal mean winds are evaluated from the so-called gradient wind expression: 


— tan0 + 2Qsin0i7 + — (2) 
a a d<p 


Values at the equator are evaluated from the curvature expression: 

1 <? 2 d> _ 

u + — ~- = l) 

2Qa d(j , 2 


(3) 


Fleming and Chandra (1989) have compared equatorial winds den ed from these data via 
Eq. 3 with monthly mean rawinsonde and rocketsonde measurements. They show that 
the derived winds reproduce the general characteristics of the quasi-biennial oscillation 
(QBO) and semi-annual oscillation (SAO), although the amplitude of the QBO is 
substantially underestimated in the derived winds. 

Zonal Fourier components of horizontal winds are derived from the zonal and 
meridional momentum equations linearized about the zonal mean zonal wind u: 
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' a "' + /v' + -L-^:=o 


a cos0 dk 
3 dv' 


a cos<p dk 


acosip dk + a~d<p 


(4a) 

(4b) 


In Eqns. (2)-(4), X is longitude, 0 latitude, u and v the zonal and meridional velocities, 
respectively, and a is earth radius (6.37 x 10^ m). Overbars denote zonal means and 
primes deviations therefrom, and 



2fl sin 6 ~(u cos <p) 

a cos0 dip J 

(4c) 

2Qsin0 +— tan^j 

(4d) 


where ft is earth rotation rate (7.3 x lO^ sec‘1). The coupled equations (4a-b) are solved 
for each zonal wave number to evaluate the spectral coefficients of u' and v'; note that 
these equations become singular for zonal wave k as 


c _ (u ■ k / acosiPr . , 

O 7TZ — 7 I 

ff 


(4e) 


For the data presented here, geostrophic winds were substituted if S > 0.5 (the results are 
not sensitive to this exact value). The use of these balanced winds and their improved 
treatment of stratospheric data (as compared to geostrophic winds) is discussed 
extensively in Randel (1987). 


Statistics of several derived quantities are also presented here. A measure of the 
zonal mean static stability is given by the Brunt- Vaisala (or bouyancy) frequency 
squared: 


N 


2 = * 


H 



( 5 ) 
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where H is the atmospheric scale height, chosen here to be 7000 m, z = H in (1000 mb/p) 
is the vertical coordinate, and K =2/7. Relatively large values of correspond to strong 
stability with respect to vertical motions. A central quantity in the quasi-geostrophic 
theory of Rossby wave propagation and stability of the zonal mean flow is the latitudinal 
gradient of zonal mean potential vorticity: 


_ 

q-y - — cos< 


1 d_ _1 d_ 

~a^ d<P cos<p d<p 

d f 


(u cosip) 


-(2Qsin#)^ e z ^ 


dz 


N 2 


( 6 ) 


A zonal cross section of potential temperature (0) and Rossby-Ertel potential vorticity 
(PV) is included for each month. Potential temperature is calculated as: 


0 = 7 


( 1000 mb 
l P 




(7) 


and potential vorticity from: 

PV = -g( 2flsin0 + O~ (8a) 

dp 

Here C is the relative vorticity, and g = 9.81 m • sec '2. For the zonal mean, this reduces 
to 

PV = -g[ 2C2sin0 “r^T(“cos0)]“7 (8b) 

a cos# d<p dp 

PV is scaled in units of 10'^ K m^ termed a PV unit Statistics of the zonally- 

averaged poleward eddy heat flux ( vT') and poleward eddy momentum flux (mV) are 
calculated and displayed. Also calculated are estimates of the zonal mean Eliassen-Palm 
(EP) flux vectors and their divergence. The orientation of the EP flux vectors indicates 
the wave group velocity or direction of wave activity flux (where such quantities are well 
defined), and the relative importance of heat versus momentum fluxes. The EP flux 
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divergence (also called the wave driving) gives a concice measure of the net effect of 
wave-induced heat and momentum fluxes in terms of a force per unit mass on the zonal 
mean flow; it is also equivalent to the poleward flux of quasi-geostrophic potential 
vorticity. A review of the motivation and theory behind EP flux diagnostics can be founu 
in Edmon et at. (1980). Here the EP flux vectors and divergences are estimated using the 
primitive equation expressions, neglecting terms involving the vertical velocity. The 
zonal mean wave driving (DF) is givei by: 


with 


DF - 


_1 

p(z) -fi COS 0 

1 


(VF) 


p(z)a cos0 


33 * £('>“■♦)♦$] 


„ , . . —j—, R du vT 

F* = p(z) • acos 0 -« V - “7 • — *- 

r [ H dz n* 


F z = p(z)-acos0 


r f i vT'i 

V « w\ 


(9) 


(10a) 

(10b) 


Here p(z) is the basic state density (proportional to e' 2 ^). EP flux diagrams here follow 
the convention in Randel et al. (1987); notably, the vectors are multiplied by e^H so 
that they are visible throughout the stratosphere. Because the use of balanced winds in 
the tropics is problematic, wave fluxes and EP diagnostics are displayed only polewards 
of 20° in each hemisphere. 


Vertical derivatives are approximated by finite differences between adjacent 
pressure levels, followed by linear interpolation in log pressure back to the original 
pressure levels. Because of the separate analysis schemes at 100 and 70 mb, and the high 
number of derivatives used, eddy heat fluxes occasionally show discontinuous behavior 
between the 100 and 70 mb levels Further vertical differentiation of these fluxes can 
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result in obviously bad EP flux divergence estimates; this problem is alleviated here by 
smoothing the eddy heat flux value? across 100 and 70 mb prior to the EP flux 
calculations. Meridional derivatives are approximated by centered finite differences on 
the Gaussian latitude grid, and longitudinal derivatives are calculated spectrally. 

2. STATISTICAL CALCULATIONS 


Zonal mean temperature and wind and zonal mean wave flux statistics presented 
here are calculated as follows ( \ffj represents the value of variable y for year i (i=l-12 

for 1979-1990) and day j (j=l-31 for January, etc.)): 

(1) Monthly means: 

M'tjSv'/ (11) 

M 


(2) Ensemble monthly means: 



( 12 ) 


(3) Ensemble average daily standard deviation (about the respective monthly means) 


a daily ~ 


12 , 31 


-.1/2 


£i£iw-w,y 

7=1 


i=V 


(13) 


(4) Interannual standard deviation of monthly means (about the ensemble average): 


^int - 


1 12 . 
TrS([^l|-<V'>) 2 

i=l 


nl/2 


(14) 


The daily and interannual variability estimates are not unbiased estimates of the variance 
(since the normalization factor is the number of samples), but have the advantage that the 
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variances Odaily and (Tint may be summed to give the total variance about the ensemble 
mean. 


Zonal averages of daily transient wave geopotential height and temperature 
fluctuations are calculated from: 


j 2n\ 

a trans = J 


ii?i 31 


il/2 


0 


i=l ;=l 


dX 


(15) 


where (ftj is the daily geopotential height or temperature and [0]/ is the monthly mean 

at each latitude, longitude and pressure. Horizontal cross sections of the daily 
geopotential height variances are also included for several pressure levels. The zonally 
averaged variability defined in Eq. IS can be compared with the standard deviation about 
the zonal mean due to all stationary waves, calc ulated from: 


G stat - 


\_ 

2 


I«s*> 2 +<c*> 2 ) 

k=\ 


(16) 


Also calculated is the ensemble average amplitude of stationary zonal waves 1-3 
individually, calculated from the twelve year ensemble average sine <Sfc> and cosine 
<Ck> coeffecients via: 

(AM/> /t ) = [(S jt > 2 +(C jt ) 2 ] 1/2 (17) 

2 

Note that the squared amplitudes (<AMPk>2 in Eq. 17) sum to twice the variance 

in Eq. 16). Stationary wave heat and momentum fluxes and EP flux divergences are 
similarly calculated from the ensemble sine and cosine coeffecients, and transient waves 
are defined as the total minus stationary contribution. Note that as defined here, 
stationary means the wave components which are constant over the entire twelve years of 
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data (and are thus in some sense truly stationary). Interannual variations in monthly 
mean wave quantites is then by definition part of the transient wave component. 

Latitude-time and height-time sections of zonal mean winds, temperatures and 
planetary wave amplitudes are calculated from three-day non-overlapping means to save 
computing space. This results in 122 values (triads) per year (except for leap years, the 
last triad is the average of days 364 and 365). These three-day mean values for the 1th 

t 

triad 0=M22) are denoted l// { - . 'Low-frequency' values are defined as 1 1 -point running 
means: 

^'=nl vi*" as) 

n=-5 


i.e., 33-day running mean (this is done to approximate monthly means). Ensemble means 
are twelve year averages of the low-frequency values: 



(19) 


Ensemble ‘daily’ variances are calculated as squared differences between the three-day 
mean and 33-day mean values: 



( 20 ) 


Interannual variances of the 33-day means is calculated from: 



( 21 ) 
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Ensemble rms zonal wave 1-3 (and 4-8) geopotential heiglu and temperature amplitudes 
are calculated from: 


rms 


■ si i K +zc u)l 

L i=i *=i v J \ 


( 22 ) 


Ensemble stationary zonal wave 1, 2 and 3 geopotential height and temperature 
amplitudes are calculated from 'low-frequency' ensemble means (eq. 18) as: 

nl/2 


'stat k 


= Hf + K} 2 


Ensemble rms transient wave amplitudes are given by: 

4ans k - jjl ,( z ^-{zs/}) 2 +(zc/ t *>{zC i f}) 2 ] 


(23) 


(24) 


with definitions as above. 

The use of three-day non-overlapping means will remove much of the variance 
with periods less than five days; this should have little effect in the stratosphere, where 
much of the variance is on time scales of 5 days and longer. Figure 1 illustrates 
December zonal mean wind daily variability calculated from daily (top) and 3-day mean 
data (bottom). Near identical patterns demonstrate the small bias introduced by using 
three day means. The use of three-day means is probably a poorer approximation for the 
RMS zonal wave 4-8 variance, because shorter scale waves tend to have faster time 
scales. 


It will be noted that there can be a significant reduction in the daily variability 
calculated from 33-day running means versus that calculated as deviations from the 
discrete monthly means (as presented in the meridional cross sections). This is due to 
excessive variance associated with using monthly means as the basis for deviation during 
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times of strong seasonal variability; an example is shown in Fig. 2. Where differences 
exist, the running mean values of triad deviations are more representative of true daily 
variability about the seasonal cycle. 'Daily variability* estimates in the seasonal diagrams 
have a spectral response function which maximizes for periods of S-SS days. 

Latitude-time sections of wave fluxes also use 122-triad values per year, but the 
fluxes are calculated for each of the individual three days prior to averaging. 'Stationary 
wave' fluxes are calculated from 33-day running mean sine and cosine spectral 
coeffecients, and 'transient wave' fluxes are simply the total minus stationary wave 
contribution. 
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Fig. 1. Meridional cros. section of the daily standard deviation of zonal mean wind in 
December (m/s), calculated using daily data (top), and 3-day non-overlapping means 
(bottom). The good agreement shows that 3-day means are appropriate for estimating 
daily variability. 




15 



33DAV 

Rl'NNINc; 

MEANS 


DAY MEANS 


Fig. 2. Time series of 3-day mean zonal wind at 1 mb, 60 S during September-November 
1979, along with the 'low-frequency' (33-day running means) and monthly mean values. 
The numbers indicate the triad standard deviations for each month, calculated about the 
monthly means, and about the running means. Note that the standard deviation about the 
monthly means are substantially larger, due to the strong seasonal transition at this time. 
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4. RESULTS 


MONTHLY AVERAGE MERIDIONAL CROSS SECTIONS 

For each month there are six pages of figures; individual plots are set in the sequence 
described in a set of templates that follow, with contour intervals noted. Note that some 
contour levels change above and below 50 mb. 
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iii:misi’IIi:ki(' polar projections for January. April. July 

AND OCTOBER 


Quantity Unite 

Monthly average geopotential heights (m) 

Monthly average doily geopotential 

height standard deviation (m) 

Monthly average wind speed isotachs 

(u^ + v2)l/2 (m/s) 

Monthly average temperature (K) 


Note: These polar stereographic plots cover approximately 15* -pole for each hemisphere, 
with latitudes of 20*. 40*. 60* and 80* noted as dashed lines. 
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ENSEMBLE AVERAGE LATITUDE-TIME DIAGRAMS 


Quantity Units 

Zonal mean temperature (K) 

Zonal mean winds (m/s) 

Zonal mean temperature daily variability (Eq. 20) (K) 

Zonal mean wind daily variability (Eq. 20) (m/s) 

Zonal mean temperature interannual variability (Eq. 21) (K) 

Zonal mean wind interannual variability (Eq. 21) (m/s) 

Meridional vorticity gradient (qy - Eq. 6) 10 -11 nr 1 s' 1 

RMS stationary and transient wave geopotential height (m) 

wave amplitudes (Eqns. 22-24) 

RMS stationary and transient temperature wave (K) 

amplitudes (Eqns. 22-24) 

Zonal wave 1-12 poleward heat flux (K-m/s) 

(total, stationary and transient waves) 

Zonal wave 1-12 poleward momentum flux (m^/s^) 

(total, stationary and transient waves) 

Zonal wave 1 - 1 2 EP flux divergence (m/s per day) 

(total, stationary and transient waves) 


NOTES: Zero contours are omitted in the wave flux diagrams. The sign of the EP flux 
divergences have been reversed for clarity (solid lines denote negative values). 
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ENSEMBLE AVERAGE HEIGHT-TIME DIAGRAMS 

Quantity UnilS 

Zonal mean temperatures (K) 

Zonal mean winds (m/s) 

Zonal mean temperature daily variability (Eq 20) (R) 

Zonal mean wind daily variability (Eq. 20) (m/s) 

Zonal mean temperature interannual variability (Eq. 21) (K) 

Zonal mean wind interannual variability (Eq. 21) (m/s) 

RMS stationary and transient wave geopotential (m) 

height wave amplitudes (Eqns. 22-24) 

RMS stationary and transient temperature (K) 

wave amplitudes (Eqns. 22-24) 


Note: The time axis has been modified for the Northern Hemisphere diagrams to run from 
July-June (i.e., NH winter is in the middle). Geopotential heights amplitude 
contours are 20 m below SO mb, and 60 m above 50 mb. 
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HEIGHT (KH1 HEIGHT I KM I HEIGHT ( KM I HEIGHT (KM ) 



PRESSURE I MB I PRESSURE rf*B l PRESSURE trei PRESSURE trj 
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LATITUDE -TIME AND HEIGHT-TIME DIAGRAMS OF ZONAL WINDS AND RMS 
GEOPOTENTIAL HEIGHT WAVE AMPLITUDES FOR INDIVIDUAL YEARS 
1979-1990 


Note: These diagrams have been constructed using three-day non -overlapping mean 
values. They are included to show seasonal cycles and actual dAily variability. 
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19/9 300 no zonal vino wy« j§» no Zonal wino i vmi j§» no zonal wind 



ja re n« At' nr in m ac sp ol nv ocja te hh ap m jn m A(. sp ut nv oc ja fit n»< ap hi jn it ac sp or m r>< 
1903 300 NB zonal WINO I MM3 JOB HU ZONAL WINO 19B« 1M Hb ZONAi WIND 



JA FB HR AP nr JN JT AC, SP OC NT OCJA FR HR AP nr IN ir A(. SP Of Nv OC JA Ffl NN AP nr JN ir AC SP OC NV IK 


)<tes ifi» nu zonal wino t<tec> 300 nu zonal wind iwbf 300 nb zonal wino 



ja fb nr ap nr in jr ac sp oc nv ocja re hr ap nr in ir ac sp oi nv oc ja fb hr ap nr jn jr ac, sp oc nv or 

I960 300 nB ZONAL WIND 1989 300 NB ZONAL WINO 1990 300 NB ZONAL WINO 



JA FB HR AP NT IN Jt AC SP Of NV OCJA I R NH AP nr JN tr AC, SP OC NV DC JA FB HR AP NT JN II AC, SP OC NV OC 





















iayitudc iayitudc latitude latitude 


220 


1R79 50 MB ZONAL WIND 1900 *0 MB ZONAL WIND 1981 be MB ZONAL WIND 



JA FB MR »P Mr JN J Y AG SP OC NV OCJA FR MR AP MY JN IT At SP OC NV DC JA FB MR AP MY JN JY AC SP OC NV DC 


W05 be MB ZONAL WIND 19Bt> be MB ZONAL UlND 1967 50 MB ZONAL WIND 



JA FB MR AP MY JN JY AC SP OC NV OCJA FB MR AP MY JN JY AC SP OC NV OCJA FB MR AP MY JN JY AC SP OC NV OC 










LATITUOt LATJYUDC LATITUOE LATtYUOt 


1979 30 MB ZONAL WIND 1980 30 MB ZONAL WIND 1981 30 MH ZONAL WIND 




JA FB NR AP NT JN JY AC SP OC NT OCJA FB HR AP MY JN JY AC SP OC NT OC JA FB MR AP NT JN JY AC SP Of NT OC 


1985 30 MB ZONAL WIND 198b 30 MB ZONAL WIND 1987 30 MB ZONAL U1NU 



1988 30 MB ZONAL WIND 1989 30 MB ZONAL WINO 1990 30 MB ZONAL WIND 



JA FB MR AP MY JN JY AC SP OC NT OCJA FB MR AP MY JN J» «C SP OC NT OCJA FB MR AP MY JN JY AC SP OC NT OC 
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t0 hu zonal v i no ivb0 10 hb zonal wind irbi 10 hb zonal vino 



jA FB HR AP nr JN JT AO SP OC NV DC JA FB HR AP HY JN JY AG SP OC NV DCJA FR HR AP HY JN JY AC SP OC NV DC 


198? IB HB ZONAL WIND 19B3 10 HB ZONAL WIND 1984 10 HB ZONAL WIND 



JA FB HR AP HY JN JY AG SP OC NV DC J* ffl HR AP HY JN JY AG SP OC NV DCjA FB HR AP HY JN JY AG SP Of NV DC 


1 98S 10 HB ZONAL WINO 1986 10 HB ZONAL UIND 1987 10 HB ZONAL WIND 



JA FB HR AP HY JN JY AG SP OC NV DCJA FB HR AP HY JN JY AG SP OC NV DCJA FB HR AP HY JN JY AG SP OC NV OC 


1 9B8 10 HB ZONAL UlNO 1909 10 HB ZONAL WINO 1990 10 HB ZONAL WIND 



JA FB HR AP HY JN JY AG SP OC NV DCJA FB HR AP HY JN JY AG SP OC NV DCJA FB HR AP HY JN JY AG SP OC NV OC 



















j» re mr ap «r jn jy ac. sp oc nv dcja fb mr ap my jn jy ag sp oc nv ocja re mr ap my jn jy ac sp oc nv 


1983 b MB ZONAL WIND 1984 b MB ZONAL W1NO 

M' ' It ’ ' ' ' I I I 


51 II 






J* re MR AP MY JN JY AC SP OC NV OCJA FB MR AP MY JN JY AC SP OC NV DCJA FB MR AP MY JN JY AC SP OC NV DC 
, QQt * hr 70NAL WIND 1986 5 MB ZONAL WIND 1,07 * " B * 0N * L MIND 
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19/9 I NH ZONAL W | NO rum I NO ZONAL yiNO 1 'IH 1 1 Hit ZONAL U I NO 



0* FB NR »P NT JN JT AC SP OC NV GCJA FB NR AP NT JN JT AC SP OC NV DC JA FB NR AP NT JN JT AG SP OC NV OC 


1 982 1 N0 ZONAL WIND 1983 1 N8 ZONAL WIND 1984 1 NB ZONAL WIND 



JA FB NR AP NT JN JT AC SP OC NV OCJA FB NR AP NT JN JT AC SP OC NV OCJA FB NR AP NT JN JT AC SP OC NV OC 


I960 1 NB ZONAL WIND 1<»B9 1 NB ZONAL WIND 1990 1 NB ZONAL WIND 



JA FB NR AP NT JN JT AC SP OC NV OCJA FB NR AP NT JN JT AG SP OC NV OCJA FB NR AP NT JN JT AG SP OC NV DC 

























ja fb mh ap my jn iy ag sn nc nv nrjA fr mr ap my in iy Ar, sp oc nv nr ja fr mr ap my jn iy ag sp nc nv nr 


19B2 I MB RMS WAVE AMP 1983 1 MB HNS WAVl AMP 1984 1 MB RMS WAVl AMP 



JA FB MR AP MY JN JY AG SP OC NV 0CJA FB MR AP MY JN JY AG SP OC NV OCJA FB MR AP MY JN JY AG SP OC NV DC 


19BS 1 MB RMS WAVC AMP 19Bt> 1 MB RMS WAVL AMP 1987 1 MB RMS WAVL AMP 



JA FB MR AP MY JN JY AG SP OC NV OCJA FB MR AP MY JN JY AG SP OC NV OCJA FB MR AP MY JN JY AG SP OC NV OC 


1988 1 MB RMS WAVL AMP 1989 1 MB RMS WAVL AMP 1990 1 MB RMS WAVE AMP 



JA FB MR AP MY JN JY AG SP OC NV DC JA FB MR AP MY JN JY AG SP OC NV OCJA fB MR AP MY JN JY AG SP 0C NV DC 
















JA FB HR AP MT JN JY AC SP OC NV OC JA FB HR AP HY JN JY AG SP OC NV DC JA FB HR AP HY JN JY AC SP OC NV OC 









latitude l»tituoe latitude 




1RB5 EQTR ZONAL WIND 1906 EQTR ZONAL WIND 1QB7 EQTR ZONAL WIND 
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TABLES OF MONTHLY MEAN ZONAL MEAN TEMPERATURES AND WINDS 
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